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Abstract 

A theory of non-local linear ion-temperature-gradient (ITG) drift modes while retaining non- 
adiabatic electrons is presented, extending the previous work [S. Moradi, et al Phys. Plasmas 18, 
062106 (2011)]. A dispersion relation is derived to quantify the effects of the fractional velocity 
operator in the Fokker-Planck equation modified by temperature gradients and non-adiabatic elec- 
trons on the real frequency and growth rate. Solving the dispersion relation, it is shown here that 
as the plasma becomes more turbulent, it deviates from a Maxwellian distribution and becomes 
Levy distributed. The resulting Levy distribution of the plasma may thus significantly alter the 
transport. The relative effect of the fractional derivative is larger on the real frequency than on 
the growth rate of the ITG mode. 

PACS numbers: 52.25.Dg, 52.30.Gz, 52.35.Kt 
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INTRODUCTION 



The high level of anomalous transport in magnetically confined fusion plasmas is still an 
unresolved issue in the quest for controlled fusion. Furthermore, a deterministic description 
of intermittent events in plasma turbulence is improper due to the stochastic nature of 
the transport exhibiting non-local interactions as well as non-Gaussian probability density 
functions (PDFs). The PDFs of heat and particle flux display uni-modal non-Gaussian 
features which is the signature of intermittent turbulence with patchy spatial structure that 
is bursty in time. The turbulent behavior in magnetically confined plasmas is the main 
ingredient in the anomalously high transport of heat, particles and momentum visible in 
present days large experiments. One crucial component of the turbulent transport is the 
so-called ion-temperature-gradient (ITG) driven turbulence. The ITG turbulence is found 
to be bursty in nature where a significant part of the transport is carried by large avalanche- 
like events. More specifically, exponential scalings are often observed in the PDF tails in 
magnetic confinement experiments, and intermittency at the edge strongly influences the 
overall global particle and heat transport. In particular it may for instance influence the 
threshold for the high confinement mode (H-mode) in tokamak experiments [1]. In view 
of these experimental results, theories built on average transport coefficients and Gaussian 
statistics fall short in predicting vital transport processes. There is a considerable amount o 



experimental evidence [2|-|5| and recent numerical gyrokinetic js-9] and fluid simulations 
that plasma turbulence in tokamaks is highly non-local. A satisfactorily understandin g of 
the non-local signatures as well as the ever-present non-Gaussian PDFs of transport [lll-ll3l] 
found in experiments and numerical simulations is still lacking. 

An attractive candidate for explaining the non-local features of ITG turbulence is by 
inclusion of a fractional velocity operator in the Fokker-Planck (FP) equation 14| yielding 
a non-local description that have non-Gaussian PDFs of heat and particle flux. The frac- 
tional operator introduces an inherently non-local description with strongly non-Maxwellian 
features of the distribution function resulting in significant modification of the transport 
process. The non-locality is introduced through the integral description of the fractional 



derivative 



15 



16]. There are a number of other phenomenological studies of the effects 



of fractional derivative models. Using fractional genera 



kinetic descriptions have been developed previously 
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izations of the Liouville equation, 

n 

, 118|]. It has been shown that the 



chaotic dynamics can be described by using the FP equation with coordinate fractional 
derivatives as a possible tool for the description of anomalous diffusion [ji]] and much work 
has been devoted on investigation of the Langevin equation with Levy white noise, see e.g. 



Refs. 
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211 ]. or related fractional FP equation. Furthermore, fractional derivatives have 



been introduced into the FP framework in a similar manner 
however a study on ITG modes is still lacking. 
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231 ] as the present work 



In this paper we introduce the Levy statistics into a Langevin equation which yields a 
fractional FP description. In order to calculate an equilibrium PDF we use a model based 
on the motion of a charged Levy particle in a constant external magnetic field obeying 
non-Gaussian, Levy statistics. This assumption is the natural generalization of the classical 
example of the motion of a charged Brownian particle with the usual Gaussian statistics. 
The fractional derivative is represented with the Fourier transform containing a fractional 

n 

exponent. Here, we extend the work presented in Ref. |14[] to include the effects of finite 
temperature gradients and non-adiabatic electrons leading to a fractional description of the 
non-local effects in ITG turbulent transport in a gyrokinetic framework. We quantify the 
non-local effects in terms of a modified dispersion relation for linear ITG modes. We have 
considered a case with constant external magnetic field and a shear-less slab geometry. The 
characteristics of the ITG modes are fundamentally changed, i.e., the values of the growth 
rate and real frequency are significantly altered by the order of the fractional derivative 
a. However, the relative effect of the fractional derivative is larger on the real frequency 
than on the growth rate of the ITG mode. This is different from the results obtained in 
Ref. [14( where the growth rate was increased strongly as the plasma deviated from a = 2 
limit (Maxwellian). We have found that the basis of this difference is due to the different 
assumptions on the electrons: adiabatic and non-adiabatic electrons are considered in Ref. 
14j and in the present work, respectively. 



The paper is organized as follows: first we present the mathematical framework of the 
fractional FP equation (FFPE) which is used to derive a dispersion relation for the ITG 
modes while retaining the non-local interactions. In the next section, the deviations from a 
Maxwellian distribution function are investigated and the dispersion relation is solved. We 
conclude the paper with a results and discussion section. 
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FRACTIONAL FOKKER-PLANCK EQUATION 

Following Ref. [jjj], the FFPE with fractional velocity derivatives in shear- less slab ge- 
ometry in the presence of a constant external force can be written as: 

dF s dF s F dF s d . „ , ^d a F s 
ot or m s av c/v a|v| a 

where s(= e, z) represents the particle species and < a < 2. Here, the term |j^§ is the 

fractional Riesz derivative. The diffusion coefficient, -D, is related to the damping term u, 



according to a generalized Einstein relation 24 \ 



l (2 ) 

ri + akr 1 



Here, T a is a generalized temperature, and force F represents the Lorentz force (due to a 
constant magnetic field and a zero-averaged electric field) acting on the particles of species 
s with mass m s and T(l + a) is the Euler gamma function. The solution of the equation 



JT]), i.e. the generalized equilibrium distribution, for a general a can be obtained as [14j : 



where 



F ( r V) - n * (r) / ^ dk l C -*W vx+kffvii) -gflki r+|kyn f3 s 



X? = ^> — f4) 



and we have introduced a generalized thermal velocity as 



V« s = (5) 



Using the generalized equilibrium distribution expressed in equation (E]), we will now quan- 
tify the non-local effects on drift waves induced by the fractional differential operator by 
determining the dispersion relation for ITG driven drift modes. We start by formulating 
the linearized gyro-kinetic theory where the particle distribution function, averaged over 
gyro-phase is of the form (see Ref. H) 

/ s (r,v) = F s (r,v) + (27r)- 4 x J J dk duuexp(ik ■ r - iut)8f£jv). (6) 

We assume that the turbulence is purely electrostatic and neglect magnetic field fluctuations 
(<5B = 0). For small deviations from the local equilibrium we find the linearized gyro-kinetic 
equation of the form 

(d t + ik\\v\\)6f£(v\\,v ± ,t) = i[—k y V x + — fc||a||]F s (x,f||,w ± )Jo(|fi s |" 1 A; ± w ± )(50 k (t). (7) 

r> trie 



Here d» = d/dvn. Evaluating explicitly the derivatives of the distribution function in equa- 
tion ([3]), we obtain the following relations: 



^k y V x F s (x,v) = y 



s dlnn s (x) 1 d In T 3i0l (x), _, . , e s s 

2 dx W^) + tT U ^ X 



{ 



n s (x) 



where 



2ti^ 2 V 7 XDJ (2tt) 3 / 2 



d In A(x) 
dx 



I T2^ [ -? (|k - |Q + |k|r) ^ n ^ (x) ]e-^ v ^ k l v ^-^' k "'" + ' k l'")}, (8) 



i d M*) anr | 

A(x) dx ' an(1 



-knduFJx, v) 



n s (x) 



m. 



m. 



l 2ir 3 / 2 VWJ (2tt) 3 / 2 



dk v t dkr. 



(9) 



where a;* k = ^gk y , and we assumed that the space dependence of F s is only in the x direction 
perpendicular to the magnetic field as well as for the density gradient. In the equation above, 
Jo is the Bessel function of order zero, v\\ is the parallel velocity, v± = (v% + v 2 ) 1 ^ 2 is the 



perpendicular velocity and hence we write the total speed as v 



+ U 2 ) 1 / 2 . The linearized 



gyro-kinetic equation could be further Laplace transformed. The Fourier-Laplace transform 
of the fluctuating electrostatic potential is 

dte^Sfait). (10) 

Similar formula defines the Fourier-Laplace transform of 5/k,w Therefore the Fourier- 
Laplace transformed gyro-kinetic equation ([7]) is 

— — fc M x7|,)^y^(v„, t7_l, *) = — V-l)*^.^ -H <5-^C«|| a 0)- (H) 

Its solution is 

^k,J^||,^)=^k,J^lh^){-Ak,J^||,^)¥k^ + ^h,^,0)}, (12) 
where the operator 

&>||,«±) = 



1 



, , -T (13) 
— i(u — k\\v\\) 

is the unperturbed propagator of the gyro-kinetic equation, and we have introduced the 
function A£. w (v\\, v±) as 

. e„ dlnn s (x) 1 d In T Si<x {x) , x . 

[ j o j \ F s\X, V)j (|"a| k ± Vs_)+l- f 



—l- 



-UJ 



*kL 



dx 



f n 8 (x) f dkldkl 
Utt^Vzd' (2tt) 3 / 2 
n s {x) 



dx 

^Vifott | a _|_ Ik" I") — — r ^ s,a ^ x ^ 
a ± " da; 



— /ci 



2^/ 2 ^XDJ (2tt) 3 / 2 



•fk^e" <(k:i - v - L+k ii V||) e"- (|k3 - |a+|k I l " ) 



i ( k I v -+ k l v ll) e -?(l k Il a +l k ll")|j (|fi s |-i^) 
Jod^r 1 ^). (14) 



Here, the wave vector perpendicular to magnetic field is k± = [k 2 x + fc^) 1//2 . The gyro- kinetic 
Equation fl6]) is complemented with Poisson equation for the electric potential. For fluctu- 
ations with wave vectors much smaller than the Debye wave vector, the Poisson equation 
becomes the quasi-neutrality condition 

E e ^ = 0, (15) 

s 

where the density fluctuation is related to the distribution function through 

s K,cu = -Y n ' S< hv> + / dvJ o{\^s\~ 1 k±v±)5f^(vi\,v ± ). (16) 

In the above equation we have separated the adiabatic response (first term on the right hand 
side) from the non-adiabatic response (second term on the right hand side). We have to 
keep in mind that the density n s coming from the F s (x, v) in the adiabatic response is also 
given by Equation ([3]) and for a general < a < 2 the adiabatic response can be different 
than that calculated by Maxwellian distribution. Using the quasi-neutrality condition (Q 
we find the dispersion equation which determines the eigenfrequencies as a function of the 
wave vector, u> = uj(k) = u r (k) + «7(k). In the simplest case we consider a plasma consisting 
of electrons and a single species of singly charged ions with equal temperatures. For the 
density fluctuation therefore we have 

K« = -n s {x)^8<j>^[M ad ' s + M£J. (17) 



Therefore, the dispersion equation as in the Ref. [25j is 

M ad,e + M e^ = _ M <id,i _ M * ^ ^g) 

where 

M ad,s = [ dv _ / ^^e- 4 ^^-^^^-^^- 1 ^ 1 ^ 1 "^^) 

J 27r^(T(l + a))- 1 /\fZVf s J (2tt) 3 / 2 y J 



gives the adiabatic contribution, and 

MZ,, = ^—J rfv^(^|,^)A^( V ||,^)J (|^r 1 ^ ± ; 



^ s dlnn s (x) l dlnT s<a (x) ^ r ^ J$(\Q, S \ 1 k ± v±) ^ I 



dx 2 dx 'J w-h^vw 2tt 3 / 2 (T(1 + a))-y\/2V 1 a 



dkldk| 

(2tt) 3 / 2 



J LO — K\\V\\ 



a 

T.s 



w-Vll 27r 3 / 2 (r(l + a))- 1 /2y2^ 

^rfkf Vt,s n k v r Q + | ^ la) djn7^) i(kl Vx+k | V|| )g _ _^£_ ( | k « + 

(27r) 3 / 2 T(l + a)a " dx 

T s , Q _ k r dv J2{\n.\- 1 k ± v ± ) 1 



u-k\\v\\ L 27r 3 / 2 (r(l + a))- 1 / 2 ^W 1 



a 

T.s 



dk v , dkf. 



Vr 



(2tt) 3 / 2 

(20) 

gives the non-adiabatic contribution. 

The analytical solutions for integrals over with an arbitrary a in the Equations (fl9l) 
and (T2"U|) requires rather tedious calculations. Instead we consider an infinitesimal deviation 
of the form a = 2 — e, where < e C 2 and expand the terms depending on a in the 
Equation (fl9|) around e = as follows 

1 =P -r<&(i k ir+i k r) = 



27r 3 / 2 (r(l + a))- 1 /2 2Vf 



T.s 



e -^(lkll 2 +lkll 2 ) 



where 



+ A(kl,k|f)e + 0[e 2 ], (21) 



^(|k-| 2 +K| 2 ) 



HK, k{) = e 8^75^7 + 27s + 21og ^ - 2 ^,J|kl| 2 + |k|| 2 ] + 

7^ T 2 iS (|kH 2 + |k|| 2 ) + V2,.(|ki| 2 tog[|kl| a ] + |k|j| 2 log[|m 2 ]) + ^ T 2 >g[F TiS ](|kl| 2 + |k||| 2 )} 

(22) 

and in Equation (1201) the expansion for the second term on the RHS gives 



v. 



T ' s f \k v { \ a + \K\ a ) V T.s 



2vr 3 / 2 (r(l + a))- x /*J2V£ a 
e-^ (|kl|a+|k J |a V Tfa (|kI| 2 + |k»| a ) 



in 



3/2 



+ £(kl,kjQe + (9[e 2 ], (23) 



where 



mi, kj) = ^-rr 2 ^{5(|kl| 2 + |kf | 2 ) - (2 lE + 21og[V T , s ])(|kl| 2 + |kf | 2 ) 
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-4(|ki| 2 log[|ki|] + |ki| 2 log[|k||]) + (-2V T>S + 7^ + log^]^ 2 J(|kl| 4 + |kjf| 4 ) 
-AV T , S + 2 lE V" + 2log[V TtS }Vl s )(|kl| 2 |m 2 ) + VS. (|kl| 4 log[|kl|] + im 4 log[|k||]) + 



^(|kll1k|| 2 )(log[|kl|] + log[|k{ll)X24) 

Here, we have used the Euler-Mascheroni constant je ~ 0.57721. 

Inserting the zeroth order terms in e from the expansion (}2~Tj) into Equation ( TT9|) will 
produce the Maxwellian adiabatic response 

M ad ' s = 1, (25) 

and by inserting the zeroth order terms in e from the expansion ff22l into Equation (120]) will 
produce the Maxwellian non-adiabatic response 



where 



= -„ / dv\\ / cfo^ ^ (I^I'V^e ^ , (26) 

V7rVy s J-00 Jo — oj + || 



T dlnn s (x) , v \\+v± 3 dlnT s (x) 

^ h,^) = -* k [^^ + - I" (27) 

By using the expansion defined by the expressions f l2~Tj) and ( 1231) to first order in e from 
Equations (TT9l) and ( 1201) . the adiabatic and non-adiabatic parts of the dispersion relation 

M ad,s and 

follows 

/oo /*oo 
dull / dv ± v ± x 
-00 jo 

/ T2^ e " (klv±+kiv " )A(k i' k ^ ))e = 1 + eWad,s - 



(2tt) 3 /2 

(28) 



and 



2 1 2 



^-i || -L 

= ^ f dt,,, r^ / B ^^ ll '^ ) J 2 (iO s l-W)e"^+e 

V^VTs "' — w + ^ll^ll 

d/nn s (x) l(ilnT s (i) P /•«> Jod^l" 1 ^^) 

{ - 2 ™* k[ ^ 2—dT- ] L^J dv±v± ~ ^ 



u-k\\v\\ 

l^ e - i(tl v i+tiV|l , A(klikI) + 



(2tt) 3 /2 
° d / 

da; 7-00 Jo ^ u — k\\v\\ 



. s dlnT s (x) r°° too JZdn^kxV^ 

27rw *k "H^ / ^1! / tJVL«L— ^ X 



2n — feii / rfvn / dv±v±-^ — — x 

m s J-oo Jo uj — k\\v || 



DISPERSION EQUATION 



We will now turn our attention to the problem of solving the dispersion relation described 
by Equation (fT8l) . In order to solve this dispersion equation we use the method described in 
Ref. [14J, where the dispersion relation is in the form 



(1 + K,J + te{W ad > e + WU = -(1 + A^J - u(W^ + WIJ. (30) 

Note that we have expanded in e e and e« for electrons and ions, respectively and that there 
exist a relation between the two see Ref. [ijj]. The first terms on the right and left hand 



sides generate the usual contributions to the dispersion equation as in Ref. [25( and the 
terms proportional to e generate the non-Maxwellian contributions. For the non-adiabatic 
Maxwellian response we have 

Ni u = 4 s Tdw f^[ W -^ r(M ' w) ]jg(Me-(^ a ), (31) 

•y 7T J-oo JO ' W — UJ 

with 

uf k (u, w) = u4[l + (u 2 + w 2 - ^) Vs }. (32) 

Here, b s = fe±Vr, s /^s, {w, u] = {v \\/Vt jS , v±/Vt, s }, we have introduced the following notation 
La = d ln d x^ ; Vs = Lx/L n and = f^k y /L n . Bar denotes normalization to k||Vr jS . The 
effects of the fractional velocity derivative can result in the non-Maxwellian contribution of 
the form 



Wl„ = ^r dw r duu[ W 1 V> W) ~ "3 ^ W V g(Me- (W) , (33) 
y7T J-oo jo w — to 



where 



tttliu, w) = u4[l - -VsMu, w) - w^tf (u, w). (34) 
The functions $(u,w), ^f(u,w) and T («,«;) are given in Appendix A. 



RESULTS AND DISCUSSION 



In this section we present the solutions of the dispersion Equation ( 1301 using Equations 

f[3"Tj) and 032]). We can find an expression for e e as: 

2 + + Nj^ 
6e 67.32 + ^^+1.42^' 1 ' 

Here, we have used the results shown in Ref. [14j: e T - — 1.42e e , W* d = 33.724, and W" d = 

23.6591, V T , e = 5.93 10 9 [cm/s], V T ,i = 1-38 10 8 [cm/s], and 6j = 0.42. We normalize all 

the frequencies to |fe|||Vr e an d we solve Equation (1351) for given values of 7 and u where 

uj = u r + ry corresponding to the real and imaginary (also called growth rate) parts of the 

eigenvalue. 

Figure [1] shows the deviation factor e e , as defined in equation (135]) calculated for given 
values of 7 and oj with m = 5 and rj e = 0. As seen in this figure, the deviation factor increases 
as the frequency and growth rate of the ITG mode increase. These results are in agreement 
with results given in Ref. jl^ . where it was shown that as the growth rate increases, e.g. 
the plasma become more turbulent, the plasma starts to deviate from a Maxwellian, and 
becomes Levy distributed. This qualitative behavior is observed for all relevant values of 
the temperature gradient through the parameter rji. Note that the relative effect on the real 
frequency is larger compared to the effect on the growth rate. This behavior is different 



from the results shown in Ref. 14] where the main effect of the deviation of the plasma 
from Maxwellian was observed on the growth rate of the density gradient mode. From our 
findings we expect that the basis of this difference is due to the difference in the assumptions 



on the electron dynamics: adiabatic or non-adiabatic electrons were assumed in Ref. 141 and 
here, respectively. However, these results in agreement with [14J suggest that as the plasma 
becomes more turbulent, it starts to deviate from a Maxwellian distribution and becomes 
Levy distributed. The resulting Levy distribution of the plasma may thus significantly 
alter the transport. Therefore, the impact of the redistribution of the plasma with different 
statistical properties has to be taken into account when calculating the transport effects. 
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FIG. 1. e e as function of to and 7, where the frequencies are normalized to |/c|||Vr e - 
APPENDIX A 

The functions $(it, w), ty(u,w) and T(u, w) are defined as followes: 

$(u,w) = ^-{wErfiNd] 

8\u\ 

((-1 + 3w 2 )(-2 lE + 21og[V T , s ]) + e w \F 1 ^[^, I -ti; 2 ]) 
-i|u|(14 - 87^ - 8m 2 + 4>y E u 2 - 8w 2 + 4>y E w 2 - 41og[V T)S ] 
-2e« 2 1 F 1 (W[^ 1 - M 2 ] - 2e™\F 1 ^>^, I, - W 2 ] 

-*Erfl[|«|]((-l + 2u; 2 )(-2 + 7s + 21og[V r J) + e w \F 1 ^[~, ~ - W 2 ]))}, (A36) 

and 

tf(it,w) = -^-r{zuErfi[|u|] 

((3 - 12^ 2 + 4^ 4 )(-8 + 3 lE + 61og[V T , s ]) - 9e w \F 1 ( 1 >°'°\~, ~ -«; 2 ]) 
+|u|(-84 + 487s + 180m 2 - 96-f E u 2 - 32m 4 + 247 S m 4 + 
180m; 2 - 96^ E w 2 - 192u 2 w 2 + 487£W 2 w; 2 - 32w 4 + 
247SW 4 + 24iEr%] - 9i-f E BrG.[u\ - 9Qiw 2 Erfi[u] + 
36z7 £ u; 2 Erfi[n] + 32zw 4 Er%] - 12ij E w 4 Brfi[u) + 241og[Vr, s ] - 
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24 M 2 log[Vr, s ] - 24u> 2 log[\/ T J - 18iErfi[u]log[y r , s ] + 72iw 2 Ev^[u]\og[V T , s ] - 
24i W 4 Er%]log[\/ T J - 6(-l + 2u 2 )(-l + 2^ 2 )( 1 F 1 ( 1 '°'°) [0, ^u 2 } + 1 F 1 ^[0, ^,w 2 ]) - 

l2e u \F^[\, I -u 2 } - l2e w \F^[\, ±, -i, 2 ] + 

gie^ErfiH^! 1 -i^ 2 ] + Bxi^p), - 24u 2 1 F 1 W)[0, ^ , it 2 ] - 

W^xF^ 1 ' - 1 ) [0, \,u 2 ]+ A%u 2 w\F^ [0, 1 « 2 ] + 6^ [0, ±w 2 }- 12uW^ [0, \ )W 2 ] - 
24w 2 1 F l ^ [0, l,w 2 } + A8u 2 w 2 1 F 1 ^ [0, J,«; 2 ] + \2u\F^ [0, \, u 2 ] 
-2Au 2 w 2 1 F 1 ^[0,\,u 2 } + 

12w 2 1 F 1 ^[0, -,w 2 } - 24A 2 1 iV 1A2) [0, -, w 2 ])}, (A37) 
and 

T(«, «;) = ^{<«Erfl[|«|]((-3 + 2n; 2 )(-2 + lE + 21og[K T> J) + e^i^ 1 ' ' 1 ' [|, ± -w 2 ]) + 
|u|(22 - 127s - 8m 2 + 47^m 2 - 8w 2 + 47^w 2 - 

6iErfi[w] + 3i7 S Erfi[M] + 4iw 2 Er%] - 2i7 B w 2 Erfi[id. - 41og[F TiS ] + 
6iErfi[«]log[V r , a ] - Aiw 2 EvQ[u]\og[V T:S ] - 2e u \F^ [| 1 it 2 ] - 

ie w \-2i + Erfi^)^ 1 ' ' 1 )^, 1 -n; 2 ])}. (A38) 

Here, \F\[a; b; z] denoting Kummer's confluent hypergeometric function and the superscripts 
represent the derivative of the hypergeometric function with respect to its parameters, for 
example: represents the derivative with respect to the first parameter, i.e. a, 

and Erflf-u] gives the imaginary error function Erf[m]/i. 
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